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Abstract 

Mining association rules is an important technique for discovering 
meaningful patterns in transaction databases. Many different measures of 
interestingness have been proposed for association rules. However, these 
measures fail to take the probabilistic properties of the mined data into 
account. We start this paper with presenting a simple probabilistic frame- 
work for transaction data which can be used to simulate transaction data 
when no associations are present. We use such data and a real- world 
database from a grocery outlet to explore the behavior of confidence and 
lift, two popular interest measures used for rule mining. The results show 
that confidence is systematically influenced by the frequency of the items 
in the left hand side of rules and that lift performs poorly to filter random 
noise in transaction data. Based on the probabilistic framework we de- 
velop two new interest measures, hyper-lift and hyper-confidence, which 
can be used to filter or order mined association rules. The new measures 
show significantly better performance than lift for applications where spu- 
rious rules are problematic. 

Keywords: Data mining, association rules, measures of interesting- 
ness, probabilistic data modeling. 

1 Introduction 

Mining association rules [3] is an important technique for discovering meaningful 
patterns in transaction databases. An association rule is a rule of the form 
X =>■ Y, where X and Y are two disjoint sets of items (itemsets). The rule means 
that if we find all items in A in a transaction it is likely that the transaction 
also contains the items in Y. 

Association rules arc selected from the set of all possible rules using measures 
of significance and interestingness. Support, the primary measure of significance, 
is denned as the fraction of transactions in the database which contain all items 
in a specific rule [3]. That is, 



supp(A =>Y) = supp(A U Y) = 



(1) 



where cxy represents the number of transactions which contain all items in X 
and Y, and m is the number of transactions in the database. 

For association rules, a minimum support threshold is used to select the most 
frequent (and hopefully important) item combinations called frequent itemsets. 
The process of finding these frequent itemsets in a large database is computa- 
tionally very expensive since it involves searching a lattice which, in the worst 
case, grows exponentially in the number of items. In the last decade, research 
has centered on solving this problem and a variety of algorithms were intro- 
duced which render search feasible by exploiting various properties of the lattice 
(see [T3] for pointers to the currently fastest algorithms). 

From the frequent itemsets all rules which satisfy a threshold on a certain 
measures of interestingness are generated. For association rules, Agrawal ct al. 
[5] suggest using a threshold on confidence, one of many proposed measures 
of interestingness. A practical problem is that with support and confidence 
often too many association rules are produced. One possible solution is to use 
additional interest measures, such as e.g. lift [5], to further filter or rank found 
rules. 

Several authors [HI HI H3 [T] constructed examples to show that in some cases 
the use of confidence and lift can be problematic. Here, we instead take a look at 
how pronounced and how important such problems are when mining association 
rules. To do this, we visually compare the behavior of support, confidence and 
lift on a transaction database from a grocery outlet with a simulated data set 
which only contain random noise. The data set is simulated using a simple 
probabilistic framework for transaction data (first presented by Hahsler et al. 
[IT] ) which is based on independent Bernoulli trials and represents a null model 
with "no structure." 

Based on the probabilistic approach used in the framework, we will de- 
velop and analyze two new measures of interestingness, hyper-lift and hyper- 
confidence. We will show how these measures are better suited to deal with 
random noise and that the measures do not suffer from the problems of confi- 
dence and lift. 

This paper is structured as follows: In Section [2] we introduce the proba- 
bilistic framework for transaction data. In Section [3] we apply the framework 
to simulate a comparable data set which is free of associations and compare the 
behavior of the measures confidence and lift on the original and the simulated 
data. Two new interest measures are developed in Section [4] and compared on 
three different data sets with lift. We conclude the paper with the main findings 
and a discussion of directions for further research. 

An implementation of the probabilistic framework and the new measures of 
interestingness proposed in this paper is included in the freely available R ex- 
tension package arules |16fM 
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Figure 1: Transactions occurring over time following a Poisson process. 

2 A simple probabilistic framework for transac- 
tion data 

A transaction database consists of a series of transactions, each transaction 
containing a subset of the available items. We consider transactions which are 
recorded during a fixed time interval of length t. In Figure [l] an example time 
interval is shown as an arrow with markings at the points in time when the 
transactions denoted by Tr\ to Tr m occur. For the model we assume that 
transactions occur randomly following a (homogeneous) Poisson process with 
parameter 9. The number of transactions m in time interval t is then Poisson 
distributed with parameter 9t where 9 is the intensity with which transactions 
occur during the observed time interval: 

P(M = m) = e ~ 6t ( 9 p m ( 2 ) 
m! 

We denote the items which occur in the database by L = {h, 1%, . . . , l n } with 
n being the number of different items. For the simple framework we assume 
that all items occur independently of each other and that for each item l± S L 
there exists a fixed probability Pi of being contained in a transaction. Each 
transaction is then the result of n independent Bernoulli trials, one for each 
item with success probabilities given by the vector p = (pi,p2, ■ ■ ■ ,p n )- Table [I] 
contains the typical representation of an example database as a binary incidence 
matrix with one column for each item. Each row labeled Tr\ to Tr m contains 
a transaction, where a 1 indicates presence and a indicates absence of the 
corresponding item in the transaction. Additionally, in Table [l] the success 
probability for each item is given in the row labeled p and the row labeled c 
contains the number of transactions each item is contained in (sum of the ones 
per column). 

Following the model, c,-, the observed number of transactions item li is con- 
tained in, can be interpreted as a realization of a random variable CV Under 
the condition of a fixed number of transactions m, this random variable has the 
following binomial distribution. 

P(d = a\M = m) = Wpf(l - Pi ) m - Ci (3) 
However, since for a given time interval the number of transactions is not 

X R is a free software environment for statistical computation, data analysis and graph- 
ics. The R software and the extension package arules are available for download from the 
Comprehensive R Archive Network (CRAN) under http://CRAN.R-project.org/ 
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Table 1: Example transaction database with transaction counts per item c and 
items success probabilities p. 



fixed, the unconditional distribution gives: 

oo 

P{Q = a)= p ( c * = *\ M = ™) ■ p ( M = m ) 

TYl — Ci 

— '"Vc, ..„-<:, e-"{et) m 



m = r>- ^ ' 



c,! ^ (m-Cj)! 

m — Ci 



The term X)m= c ^ (m-c^)! — * ™ ^ e secon d to last line in Equation J^tj is 
an exponential series with sum e' 1_Pi - )et . After substitution we see that the 
unconditional probability distribution of each Cj follows a Poisson distribution 
with parameter Pi6t. For short we will use — pidt and introduce the param- 
eter vector A = (Ai, A2, . . . , A n ) of the Poisson distributions for all items. This 
parameter vector can be calculated from the success probability vector p and 
vice versa by the linear relationship A = pOt. 

For a given database, the values of the parameter 8 and the success vectors 
p or alternatively A are unknown but can be estimated from the database. 
The best estimate for 8 from a single database is m/t. The simplest estimate 
for A is to use the observed counts c, for each item. However, this is only a 
very rough estimate which gets especially unreliable for small counts. There 
exist more sophisticated estimation approaches. For example, DuMouchel and 
Pregibon [TT] use the assumption that the parameters of the count processes 
for items in a database are distributed according to a continuous parametric 
density function. This additional information can improve estimates over using 
just the observed counts. 

Alternatively, the parameter vector p can be drawn from a parametric distri- 
bution. A suitable distribution is the Gamma distribution which is very flexible 
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and allows to fit a wide range of empirical data. A Gamma distribution together 
with the independence model introduced above is known as the Poisson-Gamma 
mixture model which results in a negative binomial distribution and has appli- 
cations in many fields |21j . In conjunction with association rules this mixture 
model was used by Hahsler [TS] to develop a model-based support constraint. 

Independence models similar to the probabilistic framework employed in this 
paper have been used for other applications. In the context of query approxi- 
mation, where the aim is to predict the results of a query without scanning the 
whole database, Pavlov et al. [24] investigated the independence model as an 
extremely parsimonious model. However, the quality of the approximation can 
be poor if the independence assumption is violated significantly by the data. 

Cadez et al. |TU] and Hollmen et al. [TI5] used the independence model to 
cluster transaction data by learning the components of a mixture of indepen- 
dence models. In the former paper the aim is to identify typical customer 
profiles from market basket data for outlier detection, customer ranking and 
visualization. The later paper focuses on approximating the joint probability 
distribution over all items by mining frequent itemsets in each component of the 
mixture model, using the maximum entropy technique to obtain local models 
for the components, and then combining the local models. 

Almost all authors use the independence model to learn something from 
the data. However, the model only uses the marginal probabilities p of the 
items and ignores all interactions. Therefore, the accuracy and usefulness of 
the independence model for such applications is drastically limited and models 
which incorporate pair-wise or even higher interactions provide better results. 
For the application in this paper, we explicitly want to generate data with 
independent items to evaluate measures of interestingness. 

3 Simulated and real-world database 

We use 1 month (t — 30 days) of real-world point-of-sale transaction data from 
a typical local grocery outlet. For convenience reasons we use categories (e.g., 
popcorn) instead of the individual brands. In the available m — 9835 transac- 
tions we found n = 169 different categories for which articles were purchased. 
This database is called "Grocery" and is freely distributed with the R extension 
package arules [16] . 

The estimated transaction intensity for Grocery is m/t = 327.5 transac- 
tions per day. To simulate comparable data using the framework, we use the 
Poisson distribution with the parameter 6t to draw the number of transactions m 
(9715 in this experiment). For simplicity we use the relative observed item fre- 
quencies as estimates for A and calculate the success probability vector p by 
X/9t. With this information we simulate the m transactions in the transaction 
database. Note, that the simulated database does not contain any associations 
(all items are independent), and thus differs from the Grocery database which 
is expected to contain associations. In the following we will use the simulated 
data set not to compare it to the real-world data set, but to show that interest 
measures used for association rules exhibit similar effects on real-world data as 
on simulated data without any associations. 

For the rest of this section we concentrate on 2-itemsets, i.e., the co-occurrences 
between two items denoted by U and lj with i,j — 1,2, ...,n and i ^ j. Al- 
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(a) simulated (b) Grocery 

Figure 2: Support distributions of all 2-itemsets (items are ordered by decreasing 
support from left to right and front to back). 



though itemsets and rules of arbitrary length can be analyzed using the frame- 
work, we restrict the analysis to 2-itemsets since interest measures for these 
associations are easily visualized using 3D plots. In these plots the x and y-axis 
each represent the items ij and lj ordered from the most frequent to the least 
frequent from left to right and front to back. On the z-axis we plot the analyzed 
measure. 

First we compare the 2-itemset support. Figure [2] shows the support distri- 
bution of all 2-itemsets. Naturally, the most frequent items also form together 
the most frequent itemsets (to the left in the front of the plots). The gen- 
eral forms of the two support distributions in the plot are very similar. The 
Grocery data set reaches higher support values with a median of 0.000203 com- 
pared to 0.000113 for the simulated data. This indicates that the Grocery data 
set contains associated items which co-occur more often than expected under 
independence. 



3.1 The interest measure confidence 

Confidence is defined by Agrawal et al. [3] as 

r (Y ^ v , supp(X U Y) 
confer) = suMx) , (5) 

where X and Y are two disjoint itemsets. Often confidence is understood as an 
estimate of the conditional probability P(Ey\Ex), were Ex (Ey) is the event 
that X (Y) occurs in a transaction [18] . 

From the 2-itemsets we generate all rules of the from ij => lj and present the 
confidence distributions in Figures [3j Confidence is generally much lower for 
the simulated data (with a median of 0.0086 to 0.0140 for the real-world data). 
Finding higher confidence values in the real-world data, which are expected to 
contain associations, indicates that the confidence measure is able to suppress 
noise. However, the plots in Figure [3] also show that confidence always increases 
with the item in the right hand side of the rule (lj) getting more frequent. This 
behavior directly follows from the way confidence is calculated. If the frequency 
of the right hand side of the rule increases, confidence will increase even if 
the items in the rule are not related (see itemset Y in Equation [5]) . For the 



G 



(a) simulated (b) Grocery 

Figure 3: Confidence distributions of all rules containing 2 items. 

Grocery data set in Figure [3^b) we see that this effect dominates the confidence 
measure. The fact that confidence clearly favors some rules makes the measure 
problematic when it comes to selecting or ranking rules. 



3.2 The interest measure lift 

Typically, rules mined using minimum support (and confidence) are filtered or 
ordered using their lift value. The measure lift (also called interest [5]) is defined 
on rules of the form X =>■ Y as 

A lift value of 1 indicates that the items are co-occurring in the database as 
expected under independence. Values greater than one indicate that the items 
are associated. For marketing applications it is generally argued that lift > 1 
indicates complementary products and lift < 1 indicates substitutes [6j [20] . 

Figure [4] show the lift values for the two data sets. The general distribution 
is again very similar. In the plots in Figures [4^a) and |4jb) we can only see 
that very infrequent items produce extremely high lift values. These values are 
artifacts occurring when two very rare items co-occur once together by chance. 
Such artifacts are usually avoided in association rule mining by using a minimum 
support on itemscts. In Figures |4jc) and (4^d) we applied a minimum support 
of 0.1%. The plots show that there exist rules with higher lift values in the 
Grocery data set than in the simulated data. However, in the simulated data 
we still find 50 rules with a lift greater than 2. This indicates that the lift 
measure performs poorly to filter random noise in transaction data especially 
if we are also interested in relatively rare items with low support. The plots in 
Figures Qc) and[4^d) also clearly show lift's tendency to produce higher values 
for rules containing less frequent items resulting in that the highest lift values 
always occur close to the boundary of the selected minimum support. We refer 
the reader to [5] for a theoretical treatment of this effect. If lift is used to rank 
discovered rules this means that there is not only a systematic tendency towards 
favoring rules with less frequent items but the rules with the highest lift will 
also always change with even small variations of the user-specified minimum 
support. 
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(a) simulated (b) Grocery 




(c) simulated with supp > 0.1% (d) Grocery with supp > 0.1% 

Figure 4: Lift distributions of all rules with two items. 



4 New measures of interest 

In the simple probabilistic model all items as well as combinations of items 
occur following independent Poisson processes. If we look at the observed co- 
occurrence counts of all pairs of two items, k and lj, in a data set with m 
transactions, we can form an n x n contingency table. Each cell can be modeled 
by a random variable C,j which, given fixed marginal counts q and cj, follows 
a hyper- geometric distribution. 

The hyper-geometric distribution arises for the so-called urn problem, where 
the urn contains w white balls and b black balls. The number of white balls 
drawn with k trials without replacement follows a hyper-geometric distribution. 
This model is applicable for counting co-occurrences for independent items 
and lj in the following way: Item lj occurs in Cj transactions, therefore, we can 
represent the database as an urn which contains Cj transactions with lj (white 
balls) and m — Cj transactions without lj (black balls). To assign item li ^ lj 
randomly to Cj transactions, we draw without replacement a transactions from 
the urn. The number of drawn transactions which we assign item lj to (and 
thus represent the co-occurrences between Z» and lj) then has a hyper-geometric 
distribution. 

It is straightforward to extend this reasoning from two items to two itemsets 
X and Y. In this case the random variable Cxy follows a hyper-geometric 
distribution with the counts of the itemsets as its parameter. Formally, the 
probability of counting exactly r transactions which contain the two independent 
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itemsets X and Y is given by 



P(Cxy = r) = 




Note that this probability is conditional to the marginal counts cx and cy . To 
simplify the notation, we will omit this condition also in the rest of the paper. 
The probability of counting more than r transactions is 

r 

P(CxY>r) = l-Y / p (CxY = i)- (8) 

i=0 

Based on this probability, we will develop the probabilistic measures hyper- 
lift and hyper-confidence in the rest of this section. Both measures quantify 
the deviation of the data from the independence model. This idea is a similar 
to the use of random data to assess the significance of found clusters in cluster 
analysis (see, e.g., [7]). 



4.1 Hyper-lift 

The expected value of a random variable C with a hyper-geometric distribution 
is 

where the parameter k represents the number of trials, w is the number of white 
balls, and b is the number of black balls. Applied to co-occurrence counts for 
the two itemsets X and Y in a transaction database this gives 

E{C X y) = — , (10) 
m 

where m is the number of transactions in the database. By using Equation [To] 
and the relationship between absolute counts and support, lift can be rewritten 
as 

lift(X => Y) = C ° nf(X ^ r) = »7>(* uy ) . = (11 ) 
V 1 supp(Y) supp(X) supp(y) E(Cxy) { ' 

For items with a relatively high occurrence frequency, using the expected 
value for lift works well. However, for relatively infrequent items, which are the 
majority in most transaction databases and very common in other domains [28j . 
using the ratio of the observed count to the expected value is problematic. For 
example, let us assume that we have the two independent itemsets X and Y, 
and both itemsets have a support of 1% in a database with 10000 transactions. 



Using Equation 10 the expected count E(Cxy) is 1. However, for the two in- 
dependent itemsets there is a P(Cxy > 1) of 0.264 (using the hyper-geometric 
distribution from Equation [8]) . Therefore there is a substantial chance that we 
will sec a lift value of 2,3 or even higher. Given the huge number of item- 
sets and rules generated by combining items (especially when also considering 
itemsets containing more than two items), this is very problematic. Using larger 
databases with more transactions reduces the problem. However, it is not always 
possible to obtain a consistent database of sufficient size. Large databases are 
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(a) simulated 



(b) Grocery 



Figure 5: Hyper-lift for rules with two items. 



usually collected over a long period of time and thus may contain outdated in- 
formation. For example, in a supermarket the articles offered may have changed 
or shopping behavior may have changed due to seasonal changes. 

To address the problem, one can quantify the deviation of the observed 
co-occurrence count cxy from the independence model by dividing it by a dif- 
ferent location parameter of the underlying hyper-geometric distribution than 
the mean which is used for lift. For hyper-lift we suggest to use the quantile 
of the distribution denoted by Qs(Cxy)- Formally, the minimal value of the 5 
quantile of the distribution of C'xy is defined by the following inequalities: 

P(Cxy < Qs(C X y)) < S and P(C XY > Qs(C'xy)) <1-S. (12) 
The resulting measure, which we call hyper- lift, is defined as 

hyper-lift 5 (X Y) = CXY (13) 
Qs(Cxy) 

In the following, we will use S = 0.99 which results in hyper-lift being more 
conservative compared to lift. The measure can be interpreted as the number 
of times the observed co-occurrence count cxy is higher than the highest count 
we expect at most 99% of the time. This means, that hyper-lift for a rule with 
independent items will exceed 1 only in 1% of the cases. 

In Figure [5] we compare the distribution of the hyper- lift values for all rules 
with two items at S = 0.99 for the simulated and the Grocery database. Fig- 
ure [5^ a) shows that the hyper-lift on the simulated data is more evenly dis- 
tributed than lift (compare to Figure |4]in Section 3.2 1. Also only for 100 of the 



n x n = 28561 rules hyper-lift exceeds 1 and no rule exceeds 2. This indicates 
that hyper-lift filters the random co-occurrences better than lift with 3718 rules 
having a lift greater than 1 and 82 rules exceed a lift of 2. However, hyper- 
lift also shows a systematic dependency on the occurrence probability of items 
leading to smaller and more volatile values for rules with less frequent items. 

On the Grocery database in Figure [5]jb) we find larger hyper- lift values of up 
to 4.286. This indicates that the Grocery database indeed contains dependen- 
cies. The highest values are observed between items with intermediate support 
(located closer to the center of the plot). Therefore, hyper-lift avoids lift's prob- 
lem of producing the highest values always only close to the minimum support 



boundary (compare Section 3.2 1 
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Table 2:2x2 contingency table for the counts of the presence (1) and absence 
(0) of the itemsets in transactions. 



Further evaluations of hyper-lift with rules including an arbitrary number of 
items will be presented in Section [4~3] 

4.2 Hyper-confidence 

Instead of looking at quantiles of the hyper-geometric distribution to form a 
lift-like measure, we can also directly calculate the probability of realizing a 
count smaller than the observed co-occurrence count cxy given the marginal 
counts cx and cy . 

c xr— 1 

P{Cxy < cxy) = ]T P{Cxy = i), (14) 

where P(Cxy = i) is calculated using Equation [7] above. A high probabil- 
ity indicates that observing cxy under independence is rather unlikely. The 
probability can be directly used as the interest measure hyper-confidence: 

hyper-confidence (X => Y) = P(C X y < c X y) (15) 

Analogously to other measures of interest, we can use a threshold 7 on hyper- 
confidence to accept only rules for which the probability to observe such a high 
co-occurrence count by chance is smaller or equal than 1 — 7. For example, if 
we use 7 = 0.99, for each accepted rule, there is only a 1% chance that the 
observed co-occurrence count arose by pure chance. Formally, using a threshold 
on hyper-confidence for the rules X —> Y (or Y =>■ X) can be interpreted as 
using a one-sided statistical test on the 2x2 contingency table depicted in 
Table [2] with the null hypothesis that X and Y are not positively related. It can 
be shown that hyper-confidence is related to the p- value of a one-sided Fisher's 
exact test. The one-sided Fisher's exact test for 2x2 contingency tables is a 
simple permutation test which evaluates the probability for realizing any table 
(see Table |2| with Cxy > cxy given fixed marginal counts [12]. The test's 
p-value is given by 

p- value = P(Cxy > cxy) (16) 



which is equal to 1 — hyper-confidence(X => Y) (see Equation 151, and gives the 
p- value of the uniformly most powerful (UMP) test for the null p < 1 (where 
p is the odds ratio) against the alternative of positive association p > 1 [22j 
pp. 58-59], provided that the p- value of a randomized test is defined as the 
lowest significance level of the test that would lead to a (complete) rejection. 

If we use a significance level of a = 0.01, we would reject the null hypothesis 
of no positive correlation if p- value < a. Using 7 as a threshold on hyper- 
confidence is equivalent to a Fisher's exact test with a = 1 — 7. 
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Note that hyper-confidence is equivalent to a special case of Fisher's exact 
test, the one-sided test on 2 x 2 contingency tables. In this case, the p- value is 
directly obtained from the hyper-geometric distribution which is computation- 
ally negligible compared to the effort of counting support and finding frequent 
itemsets. 

The idea of using a statistical test on 2 x 2 contingency tables to test for 
dependencies between itemsets was already proposed by Liu et al. 23J. The 
authors use the x 2 test which is an approximate test for the same purpose as 
Fisher's exact test in the 2-sided case. The generally accepted rule of thumb 
is that the x 2 test's approximation breaks down if the expected counts for any 
of the contingency table's cells falls below 5. For data mining applications, 
where potentially millions of tests have to be performed, it is very likely that 
many tests will suffer from this restriction. Fisher's exact test and thus hyper- 
confidence do not have this drawback. Furthermore, the x 2 test is a two-sided 
test, but for the application of mining association rules where only rules with 
positively correlated elements are of interest, a one-sided test as used here is 
much more appropriate. 

In Figures |6ja) and (b) we compare the hyper-confidence values produced 
for all rules with 2 items on the Grocery database and the corresponding simu- 
lated data set. Since the values vary strongly between and 1, we use for easier 
comparison image plots instead of the perspective plots used before. The inten- 
sity of the dots indicates the value of hyper-confidence for the rules k => lj (the 
items are again organized left to right and front to back by decreasing support). 
All dots for rules with a hyper-confidence value smaller than a set threshold of 
7 = 0.99 are removed. For the simulated data we see that the 108 rules which 
pass the hyper-confidence threshold are scattered over the whole image. For the 
Grocery database in Figure |6^b) we see that many (3732) rules pass the hyper- 
confidence threshold and that the concentration of passing rules increases with 
item support. This results from the fact that with increasing counts the test is 
better able to reject the null hypotheses. 

In Figures |7ja) and (b) we present the number of accepted rules by the 
set hyper-confidence threshold. For the simulated data the number of accepted 
rules is directly proportional to 1 — 7. This behavior directly follows from the 
properties of the data. All items are independent and therefore rules randomly 
surpass the threshold with the probability given by the threshold. For the 
Grocery data set in Figure [7^b) , we see that more rules than expected for 
random data (dashed line) surpass the threshold. At 7 = 0.99, for each of the 
n tests exists a 1% chance that the rule is accepted although it is spurious. 
Therefore, a rough estimate of the proportion of spurious rules in the set of m 
accepted rules is n(l — "f)/m. For example, for the Grocery database we have 
n = 19272 tests and for 7 = 0.99 we found m = 3732 rules. The estimated 
proportion of spurious rules in the set is therefore 5.2% which is about five 
times higher than the a of 1% used for each individual test. The reason is that 
we conduct multiple tests simultaneously to generate the set of accepted rules. 
If we are not interested in the individual test but in the probability that some 
tests will accept spurious rules, we have to adjust a. A conservative approach is 
the Bonferroni correction |26j where a corrected significance level of a* = a/ 11 
is used for each test to achieve an overall alpha value of a. The result of using 
a Bonferroni corrected 7 — 1 — a* is depicted in Figures [8ja) and (b) . For the 
simulated data set we see that after correction no spurious rule is accepted while 
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(a) simulated (b) Grocery 

Figure 6: Hyper-confidence for rules with two items and 7 > 0.99 (items are 
ordered by decreasing support from left to right and bottom to top). 

for the Grocery database still 652 rules pass the test. Since we used a corrected 
threshold of 7 — 0.99999948 (a = 5.2 • 10~ 7 ) these rules represent very strong 
associations. 

Between the measures hyper-confidence and hyper-lift exists a direct connec- 
tion. Using a threshold 7 on hyper-confidence and requiring a hyper-lift using 
the S quantile to be greater than one is equivalent for 7 = S. Formally, for an 
arbitrary rule X => Y and < 7 < 1, 

hyper-confidence (X => Y) > 7 ^ hyper-lift s (X => Y) > 1 for 7 = S. (17) 

To prove this equivalence, let us write F(c) — P(Cxy < c) for the distribu- 
tion function of Cxy, so that the hypcr-confidcncc h = h(c) equals F(c— 1). 
We note that F as well as its quantile function Q are non-decreasing and that 
for integer c in the support of Cxy, Qf(c) = c - Hence, provided that h(c) > 0, 
Qh(c) = c— 1. What we need to show is that h(c) > 7 iff c > Q 7 . If < 7 < h(c), 
it follows that Q 7 < Qh( c ) = c ~ 1> i- e - ; c > Q-y- Conversely, if c > Q 1 , it follows 
that c — 1 > Qj and thus h(c) — F(c — 1) > F(Q 7 ) > 7, completing the proof. 

Hyper-confidence as defined above only uncovers complementary effects be- 
tween items. To interpret using a threshold on hyper-confidence as a simple 
one-sided statistical test makes it also very natural to adapt the measure to 
find substitution effects, items which co-occur significantly less together than 
expected under independence. The hyper-confidence for substitutes is given by: 

cxy 

hyper-confidence sub (X =>Y) = P(C XY > c XY ) = 1 - Yl P ( C xy = i) (18) 

*=o 

Applying a threshold ^ sub can be again interpreted as using a one-sided 
test, this time for negatively related items. One could also construct a hyper- 
lift measure for substitutes using low quantiles; however, its construction is less 
straightforward . 
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hyper-confidence threshold y hyper-confidence threshold y 

(a) simulated (b) Grocery 

Figure 7: Number of accepted rules depending on the used hyper-confidencc 
threshold. 
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(c) simulated (d) Grocery 

Figure 8: Hyper-confidencc for rules with two items using a Bonferroni cor- 
rected 7. 
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(a) simulated (b) Grocery 

Figure 9: Hyper-confidence for substitutes with 7 sub > 0.99 for rules with two 
items. 



Database 


Grocery 


T10I4D100K 


Kosarak 


Type 


market basket 


artificial 


click-stream 


Transactions 


9835 


100,000 


50,000 


Avg. trans, size 


4.41 


10.10 


3.00 


Median trans, size 


3.00 


10.00 


7.98 


Distinct items 


169 


870 


18,936 



Table 3: Characteristics of the used databases. 



In Figures |9ja) and (b) we show the rules with two items which surpass a 
7 sub of 0.99 in the simulated and the Grocery database. In the simulated data 
we see that the 68 falsely found rules are regularly scattered over the lower left 
triangle. In the Grocery database, the 116 rules which contain substitutes are 
concentrated for a few items (the line clearly visible in Figure |9jb) corresponds 
to the item 'canned beer' which has a strong substitution effect for most other 
items). As for complements, it is possible to use Bonferroni correction. 

4.3 Empirical results 

To evaluate the proposed measures, we compare their ability to suppress spu- 
rious rules with the well-known lift measure. For the evaluation, we use in ad- 
dition to the Grocery database two publicly available databases. The database 
"T10I4D100K" is an artificial data set generated by the procedure described by 
Agrawal and Srikant [4 which is often used to evaluate association rule mining 
algorithms. The third database is a sample of 50,000 transactions from the 
"Kosarak" database. This database was provided by Bodon [B] and contains 
click-stream data of a Hungarian on-line news portal. As shown in Table [3] the 
three databases have very different characteristics and thus should cover a wide 
range of applications for data from different sources and with different database 
sizes and numbers of items. 
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Database 


Grocery/ sim. 


T10I4D100K/ sim. 


Kosarak/ sim. 


Min. support 


0.001 


0.001 


0.002 


Found rules 


40943/8685 


89605/9288 


53245/2530 


lift > 1 


40011/5812 


86855/5592 


51822/1365 


lift > 2 


27334/ 180 


84880/ 


42641/ 


hyper-lift gg > 1 


30083/ 196 


86463/ 150 


51151/ 23 


hyper-lift 99 > 2 


1563/ 


83176/ 


37683/ 


hyper-conf. > 0.9 


36724/1531 


86647/1286 


51282/ 240 


hyper-conf. > 0.9999 


15046/ 1 


86207/ 


51083/ 



Table 4: Number of rules exceeding a lift and hyper-lift (5 = 0.99) of 1 and 2, 
and a hyper-confidence of 0.9 and 0.9999 on the three databases and comparable 
simulated data sets. 



For each database we simulate a comparable association-free data set follow- 
ing the simple probabilistic model described above in this paper. We generate 
all rules with one item in the right hand side which satisfy a specified minimum 
support (see Table [4| . Then we compare the impact of lift and confidence with 
hyper-lift and hyper-confidence on rule selection. In Table [4] we present the num- 
ber of rules found using the preset minimum support and the number of rules 
which also have a lift greater than 1 and 2, a hyper-lift with 5 — 0.99 greater 
than 1 and 2, or a hyper-confidence greater than 0.9 and 0.999, respectively. 
From the results in the table we see that, compared to the real databases, in the 
simulated data sets only a much smaller number of rules reaches the required 
minimum support. This supports the assumption that these data sets do not 
contain associations between items while the real databases do. If we assume 
that rules found in the real databases are (at least potentially) useful associa- 
tions while we know that rules found in the simulated data sets must be spu- 
rious, we can compare the performance of lift, hyper-lift and hyper-confidence 
on the data. In Table [4] we see that there obviously exists a trade-off between 
accepting more rules in the real databases while suppressing the spurious rules 
in the simulated data sets. In terms of rules found in the real databases versus 
rules suppressed in the simulated data sets, hyper-lifto.99 > 1 lies for all three 
databases between lift > 1 and lift > 2 while hyper- lift . 99 > 2 never accepts 
spurious rules but also reduces the rules in the real databases (especially in the 
Grocery database). The same is true for hyper-confidence with a threshold of 
0.9 the number of resulting rules lying in between the results for the two lift 
thresholds and for 0.999 hyper-confidence only once (for the Grocery database) 
accepts a single rule. 

To analyze the trade-off in more detail, we proceed as follows: We vary the 
threshold for lift (a minimum lift between 1 and 3) and assess the number of 
rules accepted in the databases and the simulated data sets for each setting. 
Then we repeat the procedure with confidence (a minimum between and 1), 
with hyper-lift (a minimum hyper-lift between 1 and 3) at four settings for S (0.9, 
0.99, 0.999, 0.9999) and with hyper-confidence (a minimum threshold between 
0.5 and 0.9999). We plot the number of accepted rules in the real database 
by the number of accepted rules in the simulated data sets where the points 
for each measure (lift, hyper-confidence, and hyper-lift with the same value for 



5) are connected by a line to form a curve. The resulting plots in Figure 10 
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hyper-confidence 
lifl 

confidence 



o lift 

—I— hyper-confidence 

-x- hyperlift S = 0.9 

-o-- hyperlift S = 0.99 

-v- hyperlift S = 0.999 

— is— hyperlift 5 = 0.9999 



1000 2000 3000 4000 5000 6000 
simulated (accepted rules) 



100 200 300 400 

simulated (accepted rules) 



(a) Grocery 





■ o ■ 


lift 




hyper-confidence 




hyperlift 5 - 0.9 




hyperlift 5 - 0.99 


-v- 


hyperlift 5 = 0.999 




hyperlift 5 = 0.9999 



1000 2000 3000 4000 
simulated (accepted rules) 



100 200 300 400 

simulated (accepted rules) 



(b) T10I4D100K 





400 600 800 1000 
simulated (accepted rules) 





lift 




hyper-confidence 




hyperlift S = 0.9 




hyperlift 5 - 0.99 


-v- 


hyperlift S = 0.999 




hyperlift S = 0.9999 



50 100 150 

simulated (accepted rules) 



(c) Kosarak 

Figure 10: Comparison of number of rules accepted by different thresholds for 
lift, confidence, hyper-lift (only in the detail plots to the right) and hyper- 
confidence in the three databases and the simulated data sets. 
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arc similar in spirit to Receiver Operating Characteristic (ROC) plots used in 
machine learning |25j to compare classifiers and can be interpreted similarly. 
Curves closer to the top left corner of the plot represent better results, since they 
provide a better ratio of true positives (here potentially useful rules accepted in 
the real databases) and false positives (spurious rules accepted in the simulated 
data sets) regardless of class or cost distributions. 

Confidence performs considerably worse than the other measures and is only 
plotted in the left hand side plots. For the Kosarak database, confidence per- 
forms so badly that its curve lies even outside the plotting area. 

Over the whole range of parameter values presented in the left hand side 
plots in Figure [10] there is only little difference visible between lift and hyper- 
confidence visible. The four hyper-lift curves are very close to the hyper- 
confidence curve and are omitted from the plot for better visibility. A closer in- 
spection of the range with few spurious rules accepted in the simulated data sets 
(right hand side plots in Figure 10 1 shows that in this part hyper-confidence and 
hyper-lift clearly provides better results than lift (the new measures dominate 
lift). The performance of hyper-confidence and hyper-lift are comparable. The 
results for the Kosarak database look different than for the other two databases. 
The reason for this is that the generation process of click-stream data is very 
different from market basket data. For click-stream data the user clicks through 
a collection of Web pages. On each page the hyperlink structure confines the 
user's choices to a usually very small subset of all pages. These restrictions are 
not yet incorporated into the probabilistic framework. However, hyper-lift and 
hyper-confidence do not depend on the framework and thus will produce still 
consistent results. 

Note that in the previous evaluation, we did not know how many accepted 
rules in the real databases were spurious. However, we can speculate that if 
the new measures suppress noise better for the simulated data, it also produces 
better results in the real database and the improvement over lift is actually 



greater than can be seen in Figure 10 



Only for synthetic data sets, where we can fully control the generation pro- 
cess, we know which rules are non-spurious. We modified the generator de- 
scribed by Agrawal and Srikant [I] to report all itemsets which were used in 
generating the data set. These itemsets represent all non-spurious patterns 
contained in the data set. The default parameters for the generator to produce 
the data set T10I4D100K tend to produce easy to detect patterns since with 
the used so-called corruption level of 0.5 the 2000 patterns appear in the data 
set only slightly corrupted. We used a much higher corruption level of 0.9 which 
does not change the basic characteristics reported in Table [3] above but makes 
it considerably harder to find the non-spurious patterns. 

We generated 100 data sets with 1000 items and 100,000 transactions each, 
where we saved all patterns used for the generation. For each data set, we 
generate sets of rules which satisfy a minimum support of 0.001 and different 
thresholds for hyper-confidence, lift and confidence (we omit hyper-lift here since 
the results are very close to hyper-confidence). For each set of rules, we count 
how many accepted rules represent patterns which were used for generating the 
corresponding data set (covered positive examples, P) and how many rules are 
spurious (covered negative examples, N). To compare the performance of the 
different measures in a single plot, we average the values for P and N for each 
measure at each used threshold and plot the results (Figure 11 1. 
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Figure 11: Average PN graph for 100 data sets generated with a corruption rate 
of 0.9. 



A plot of corresponding P and N values with all points for the same measure 
connected by a line is called a PN graph in the coverage space which is similar to 
the ROC space without normalizing the X and Y-axes [13] . PN graphs can be 
interpreted similarly to ROC graphs: Points closer to the top left corner indicate 
better performance. Coverage space is used in this evaluation since, other than 
most classifiers, association rules typically only cover a small fraction of all 
examples (only rules generated from frequent itemsets generate rules) which 
makes coverage space a more natural representation than ROC space. 

Averaged PN graphs for hyper-confidence, lift, confidence and the \ 2 statis- 
tic are presented in Figure |11| Hyper-confidence dominates lift by a consider- 
ably larger margin than in the previous experiments reported in Figure |10[ b) 
above. This supports the speculation that the improvements achievable with 
hyper-confidence are also considerable for real world databases. Using a varying 
threshold on the x 2 statistic as proposed by Liu et al. [23] performs better than 
lift and provides only slightly inferior results than hyper-confidence. 

We also inspected the results for the individual data sets. While the char- 
acteristics of the data sets vary sometimes significantly (due to the way the 
patterns used in the generation process are produced; see @]), all data sets 
show similar results with hyper-confidence dominating all other measures. 
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5 Conclusion 



In this contribution we used a simple independence model (a null model with 
"no structure") to simulate a data set with comparable characteristics as a 
real- world data set from a grocery outlet. We visually compared the values of 
different measures of interestingness for all possible rules with two items. In 
the comparison we found the same problems for confidence and lift, which other 
authors already pointed out. However, these authors only argued with specially 
constructed and isolated example rules. The analysis used in this paper gives a 
better picture of how strongly these problems influence the process of selecting 
whole sets of rules. Confidence favors rules with high-support items in the 
right hand side of the rule. For databases with items with strongly varying 
support counts, this effect dominates confidence which makes it a bad measure 
for selecting or ranking rules. Lift has a strong tendency to produce the highest 
values for rules which just pass the set minimum support threshold. Selecting or 
ranking rules by lift will lead to very unstable results, since even small changes 
of the minimum support threshold will lead to very different rules being ranked 
highest. 

Motivated by these problems, two novel measures of interestingness, hyper- 
lift and hyper-confidence, are developed. Both measures quantify the deviation 
of the data from a null model which models the co-occurrence count of two 
independent itemsets in a database. Hyper-lift is similar to lift but uses instead 
of the expected value a quantile from the corresponding hyper-geometric distri- 
bution. The distribution can be very skewed and thus hyper-lift can result in 
significantly different ordering of rules than lift. Hyper-confidence is defined as 
the probability of realizing a count smaller than the observed count and from 
its setup related to a one-sided Fisher's exact test. 

The new measures do not show the problematic behavior described for con- 
fidence and lift above. Also, both measures outperform confidence, lift, and the 
\ 2 statistic on real-word data sets from different application domains as well 
as in an experiment with simulated data. This indicates that the knowledge 
of how independent itemsets co-occur can be used to construct superior mea- 
sures of interestingness which improve the quality of the rule set returned by 
the mining algorithm. 

A topic for future research is to develop more complicated independence 
models which incorporate constraints for specific application domains. For ex- 
ample, in click-stream data, the link structure restricts which pages can be 
reached from one page. Also the generation of artificial data sets which incor- 
porate models for dependencies between items is an important area of research. 
Such data sets could greatly improve the way the effectiveness of data mining 
applications is evaluated and compared. 
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